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Abstract. We obtain the exact renormalization map and plots of Lee- Yang and 
Fisher zeros distributions for Potts models on a number of hierarchical lattices: 
the diamond hierarchical lattice, a lattice we call spider web, the Sierpinski gasket 
and cylinders. Such models are only examples among the ones we can study in 
the general framework of hierarchical lattices, developed in a previous paper. 



1. Introduction 

Spin models on hierarchical lattices are a large class of exactly soluble models that 
have been first considered as approximations to regular lattices [H El [3] and then as 
examples of lattices invariant under a real-space renormalization procedure [H |5l [6l [7] . 
The renormalization group action for such models is therefore exact and the study of 
its dynamics provides some interesting results that can be useful in studying the 
renormalization group action in more general cases. In this paper we consider some 
examples of Potts models on hierarchical lattices, namely the diamond hierarchical 
lattice (section [3]), the spider web (section |4|), the Sierpinski gasket (section [5]) and 
cylinders (section [6]). Apart from the case of the spider web which (to the best of the 
author's knowledge) has not yet been subject of research, models on the other lattices 
have been extensively studied before (e.g. [HI [HI [10] , [HI [12] ) ; however, the purpose 
of this paper is to present all such models with a consistent and uniform method 
which also allows for the presence of an external magnetic field. This approach has 
been presented in a previous paper [13] and can be applied in full generality to all 
hierarchical lattices. For each model we will write the exact renormalization group 
generator and provide numerical results for the distribution of Lee- Yang and Fisher 
zeros. Such results are obtained using techniques which we explain in section [2] and 
in the appendices. We also report some observations which arise quite naturally from 
the analysis of the aforementioned models and also provide some new results. In 
particular, we observe that Lee- Yang zeros responsible for the infinite susceptibility of 
Ising model on the diamond hierarchical lattice in the paramagnetic phase are given 
by interactions that are only finitely renormalizable. Moreover, we are able to write 
the exact renormalization map associated to the Potts model on a Sierpinski gasket 
for all values of q. We refer the interested reader to [TFJ for a detailed treatment of 
Potts models on hierarchical lattices; what follows in this introduction is an attempt 
to summarize consistently all basic concepts we need in this paper. 
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Hierarchical lattices are infinite lattices in which we allow multiple-spin connections 
and that are obtained by iterating a decoration procedure on a finite lattice; this 
procedure amounts to substitute each edge of a lattice with a given block of spins and 
edges (see e.g. figures [H5l8p . In [13] we showed that, for such models, we can define 
a renormalization map that acts as a polynomial map on the complex vector space of 
Boltzmann weights ex.-p{—(3J) associated to local interactions J. We observe that a 
different normalization of the Boltzmann weights or, equivalent ly, a different choice 
of zero of energies, does not change the thermodynamics of the models. Therefore, 
we argue that the space of Boltzmann weights can be considered as a projective 
space and the renormalization map will act on such projective space as a rational 
map. In general, if we consider models with several type of interactions, then the 
renormalization map will act on the Cartesian product of several projective spaces 
(i.e a so-called multiprojective space) that we call dynamical space. We define physical 
space to be the space of Boltzmann weights associated to interactions given by pair 
interactions and (possibly) by coupling with an external magnetic field. In general, the 
physical space is a submanifold of the dynamical space which is not preserved by the 
renormalization map. This amounts to the well-known fact that the renormalization 
of pair interactions introduces new multiple-spin interactions. Hierarchical lattices are 
such that all possible multiple-spin interactions that arise from the renormalization 
process form a finite-dimensional space; in this sense we say that hierarchical lattices 
are exactly renormalizable. 



2. Numerical approaches 

In the next sections we will perform a numerical study of rational maps that generate 
the renormalization group of some examples of hierarchical lattices. We are, in fact, 
interested in finding the distribution of Lee- Yang and Fisher zeros for such models. 
Given the renormalization map of the model, one method to obtain numerically such 
distributions is to find all basins of attraction of stable fixed points of the map; the 
boundary of such regions is going to be the unstable set for the renormalization 
map (i.e. the so-called Julia set) and phase transitions of the model will appear 
for interactions belonging to such set (see [Appendix A.l for more details). A 



second approach, in some sense more straightforward, proceeds by computing an 
approximation of a real function called Green function. This function is a purely 
dynamical object and it is related to the free energy of the model the map is associated 
to; in particular we expect the two functions to have the same domain of analyticity 
(although this fact has been formally proved only for some cases). Once we obtain the 
numerical approximation to the Green function, applying the Laplacean differential 
operator yields the density of the measure supported on the Lee- Yang and Fisher zeros 
of the model (see [Appendix A.2| for details). 

As we pointed out in the introduction, the renormalization group action on Boltzmann 
weights is generated by a rational map on a multiprojective space ^ called dynamical 
space which contains all multiple-spin interactions that can be generated by the 
renormalization process. We will often consider a submanifold ^ that we call physical 
space. This submanifold is given by Boltzmann weights associated to interaction that 
are induced by pair interactions and possibly an external magnetic field. Let us define 
the pair interactions, i.e. let Jg be the energy given to two parallel neighbouring spins 
and Jd the energy associated to two neighbouring spins that are in different states. 
The Boltzmann weights associated to the corresponding energies will be denoted by 
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[zg : Zd] and belong to the one-dimensional complex projective space P^. A magnetic 
field, if present, will assign energy to one special state among the q Potts states and 
energy to all other states. The Boltzmann weights associated to the corresponding 
energies will be denoted by [h^ : h^] G P^. 

For each hierarchical lattice we can therefore define a map from P^ x P^ ^ ^ C ^ 
that gives projective coordinates to The pair interaction Boltzmann weights will 
belong to the first P-*^ factor and the magnetic field weight to the second P^ factor. 
We are now going to define standard local charts (coordinates) on each P^ factor; all 
numerical computations will be performed in one of such charts. Notice that all the 
coordinates we are going to define are just standard (inhomogeneous) charts on the 
projective line P^. 

Definition 2.1 We call standard interaction coordinates the coordinate chart ofF^ 
given by ( = z^/zd = exp {-f3 (Jg - Jd)) (for Zd 7^ 0). 

When dealing with zero-temperature phase transitions we will need to consider the 
inverse chart, l/( (for z^ 7^ 0^; we will call such chart inverse interaction coordinates. 

Note that the standard interaction coordinates could be obtained by setting = 
and considering the Boltzmann weights corresponding to such choice of zero of 
energies. In this sense we call them standard. In such coordinates, = and 
C = 00 are respectively the antiferromagnetic and ferromagnetic points, while C = 1 
is the paramagnetic point. The latter is fixed by all renormalization maps, while the 
ferromagnetic point = 00) is fixed whenever the hierarchical lattice is connected; 
the antiferromagnetic point is usually mapped to the ferromagnetic point by the RG 
map. 

Definition 2.2 We call standard field coordinates the coordinate h = hu/hu = 
exp(-^(i7^-i7^)). 

Again, the standard field coordinates can be obtained as the Boltzmann weights 
associated to the choice = 0. In the standard field coordinates, h = 1 corresponds 
to the case with zero field, h = 00 corresponds to the case of infinite field and = 
is when the privileged state costs infinite energy and it is therefore never assumed. 

3. Diamond hierarchical lattices 

Diamond hierarchical lattices (DHLs) have been the first hierarchical lattices to be 
investigated using tools from complex dynamics [8j; they are lattices on a standard 
graph, and they can be obtained by iterating the decoration procedure illustrated in 
figure [TJ Recall that we define the interaction as Jg if two neighbouring spins are in 



b 




Figure 1. Decoration generating DHL5 (left) with some iterations of the 
decoration procedure for DHL2 on a starting graph V (right). The hierarchical 
lattice DHLfo is the limit graph that we obtain by indefinitely iterating the 
decoration procedure. 
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the same state and if they are in different states. For this lattice the dynamical 
variables with no magnetic field are the Boltzmann weights [z^ : z^] = [z^ : ^d] ^ 
relative to the states of two neighbouring spins: 

z^ = exp(— ^Jg) same state 

= exp(— ;^Jd) different states 

In this case the physical space of the model without external field coincides with the 
dynamical space. The renormalization map can be easily written in the dynamical 
variables for all values of q; let <^ and ^ are the renormalized variables; then: 

{zi+{q-l)-zl)'' (la) 

= {2-z^-Zs+{q-2)-z^)'' {lb) 

For every b the map has a fixed point at [1 : 1] (paramagnetic point) and at [1 : 0] 
(ferromagnetic point). Another fixed point appears when b is odd at [— 1 : 1]. The fixed 
point [1 : 1] is always superattracting (i.e. the map has zero derivative at the fixed 
point), while the ferromagnetic fixed point is superattracting only if 6 > 1, therefore 
excluding the one-dimensional chain case. Thus, in all other cases, and for all values of 
g, we expect a phase transition at finite temperature, since basins of attraction of an 
attracting fixed point of a rational map always contain a neighbourhood of the fixed 
point. In figure [2] we show the aforementioned basins of attraction of the ferromagnetic 
and paramagnetic fixed point for various values of b and q. 



3.1. Magnetic field 

As explained in the previous paper it is possible to deal in a completely analogous 
way with an applied magnetic field; corresponding Boltzmann weights will appear as 
parameters of the renormalization map. In this case the dynamical variables are the 
Boltzmann weights [z^ \ z^ : z^ : Zg] G relative to the states of two neighbouring 
spins according to the following rules: 

z^ same state (special) 

Z| different states (one special) 

z^ same state (not special) 

z^ different states (not special) 



(2) 



Now given J^^J^^H^^ we can define the physical space as given by: 

z^ = exp{-p{Js + 2H,)) Z| = exp{-P{J^ + i^a + H^)) 
z^ = exp(-^( Js + 2H^)) z^ = exp(-^( Jd + 

The renormalization map is given in (|3fl |) - (|3^ : 

{h,-zl^{q-l)-h,-zl)' (3a) 

^ {ha - z^- z^^ hu- z^- z^^ {q - 2) ■ hu- z^ - z^)^ (36) 

Xn = {hu • ^hu-zl^{q-2)-hu- zlf (3c) 

^ = [hu ■ z'^ ^ 2 ■ hu ' z^ ■ z^^ {q - i) ■ hu ■ z'^)^ {Zd) 

The renormalization map does not preserve the physical space, i.e. the image of ([2]) 
(that is in general a submanifold of codimension 1) unless q = 2. In fact, in such 
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Figure 2. Standard interaction coordinates: Basins of attraction of the 
paramagnetic fixed point (in blue) and the ferromagnetic one (yellow) for the 
DHL for various values of b and q; rows have respectively b = 2,4,8,16, and 
columns have q = 2, 3, 4 

case the dynamical space is given by [z^j : : z^] ^ IP^ and the map given by ([2|) is 
surjective. One could in principle write the renormalization map in terms of Boltzmann 
weights associated to Jg, Jd, i^a, i^n (see e.g. [3 [10]); however, the map obtained in 
such variables is not rational (since it involves square roots) and its analysis is not as 
straightforward as it would be on the dynamical space. In any case it is convenient to 
perform the analysis in the dynamical space and then restrict to the physical space to 
obtain plots and thermodynamical quantities. 

One can compute with good approximation the Green function of the renormalization 
map and look for phase transitions in the magnetic field part of the dynamical space. 
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Looking just at the Ising case, with no surprise we find the full Lee- Yang circle 
for the ferromagnetic phase, and we obtain an anomalous plot for the supposedly 
paramagnetic phase (figure [3] shows the b = 2 case). The anomalous plot illustrates 




Figure 3. Standard field coordinates: Lee- Yang zeros of the diamond hierarchical 
lattice with 6 = 2; the center of both figures is z. Left: zeros in the real 
ferromagnetic phase; right: the anomalous zeros for the supposedly paramagnetic 
phase 

two interesting facts. The first (proved in [9j) is that zeros of the partition function do 
accumulate on the positive real axis even in the supposedly paramagnetic phase i.e. 
the system exhibits infinite susceptibility in the paramagnetic phase, which therefore 
is more appropriately called critical phase. The critical phase nevertheless exhibits 
paramagnetic behaviour (this is also proved in [9j); in fact, we report in figure [H the 
numerical data for the spontaneous magnetization. The second interesting fact to note 




Figure 4. Spontaneous magnetization for the Diamond Hierarchical Model 
(b = 2). The horizontal axis corresponds to real values of z, on the y axis we 
have spontaneous magnetization (in arbitrary units). We notice the presence of 
the three phases: antiferromagnetic, critical (paramagnetic) and ferromagnetic. 



is that points that are accumulating towards the positive real axis in the critical phase 
are not ordinary zeros of the partition function, but are preimages of the so-called 
indeterminacy set. In fact, the anomalous zeros in figure [3] on the right, can be seen 
as xs that decrease in size as they become dense, whereas the regular zeros in figure [3] 
on the left form a solid line. The indeterminacy set is the set of points on which the 
renormalization map is not defined, i.e. the points that would map to all Boltzmann 
weights equal to under the renormalization map (see [IS]). Points accumulating on 
the positive real axis in figure [3] correspond therefore to interactions that are only 
finitely renormalizable. Such points are in some sense anomalous from the points of 
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view of both dynamics and physics, and it would be quite interesting to understand 
if this connection is more than just a mere coincidence. 

4. Spider web 

The spider web lattice is obtained by iterating the decoration & shown in figure figure[5] 
infinitely many times; as the picture illustrates this lattice is based on what is called a 




Figure 5. Decoration for the spiderweb along with some iterations on a basic 
hypergraph 

3-uniform hypergraph. In this case the dynamical variables are the Boltzmann weights 
[z^ : : Z|] G relative to the states of 3 neighbouring spins: 

z^ same state 

two in the same state, third in different state 
z^ three different states 

We can consider pair interactions given by Jg, Jd on each edge of each triangle in the 
following way: 

z^ = exp{-(3 ■ 3Js) = exp(-/?(Js + 2Jd)) Z| = exp(-/3 • 3Jd) (4) 

which follows by giving to each dynamical variable the Boltzmann weight associated 
to the energy of the pair interactions in the corresponding configuration. In this 
case, each side of each triangle (apart from the three sides of the initial hypergraph) 
is counted twice, as each side is shared by two 3-edges. Since this multiplicity is 
uniform for (almost) all sides, this is not an issue; the renormalization transformation 
is therefore easily written in the dynamical variables: 

^^ = ^i + (g-l)-4 (5a) 
= z^' z^^ z^^ {q-2) ■ z^' z| (56) 

5| =3-4-^| + (g-3)-z| (5c) 

Notice that, in general, the renormalization map does not preserve the physical 
space submanifold. Once more, this amounts to the well-known fact that in general 
renormalizing pair interactions gives rise to interactions that cannot be written as pair 
interactions. This did not happen in the previous case because the DHL is naturally 
defined using only 2-edges. Notice, moreover, that if q = 2, the equation for z^ 
uncouples from the first two and we have that the projective space generated by the 
first two variables is invariant under the renormalization map. This is not unexpected 
since, if g = 2, there cannot be a configuration for which all three spins are in different 
states. As a matter of fact it is interesting to compute the restriction of the map in 
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such a case, as we recover a map of the quadratic family best known as the cauliflower 
(see for example [H]): 



C =r] = — 



cm 
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Figure 6. Standard interaction chart: Fisher zeros for q = 2,3,4; as we change 
q we can track the evolution of the antiferromagnetic phase. 



We can recognize the cauliflower in the physical variables in the leftmost picture of 
figure m From the map we easily see that the paramagnetic point [1 : 1] is a parabolic 
fixed point (i.e. its multiplier is a root of unity.) and the convergence of the Green 
function in its neighbourhood is rather slow. The same slow convergence rate can also 
be noticed for all preimages of this point . All points inside the cauliflower (therefore all 
antiferromagnetic interactions) will converge to the paramagnetic fixed point, while all 
points outside will converge to the ferromagnetic fixed point at infinity. This could be 
explained by the fact that frustration prevents the formation of an antiferromagnetic 
phase. For q = 3 the map acts on the full P^. This map has one indeterminacy point 
at [0:0:1], represented by the cross at the center of the appropriate plot in figure O 
In fact, this point corresponds to an interaction that allows only for the configuration 
given by all three spins in different states. It is easy to check (see figure [5j) that if 
g = 3 it is not possible to have a configuration of spins on the spiderweb satisfying 
this requirement. To this extent, this interaction is not renormalizable. The newborn 
region that surrounds the indeterminacy point is mapped to the basin of attraction 
of the ferromagnetic point oo, indicating that the behaviour of this phase could be 
antiferromagnetic. Considering q = 4 or higher we observe that the antiferromagnetic 
phase disappears. 

It is straightforward to write the renormalization group map in presence of an external 
magnetic field. However, for sake of clarity, we restrict ourselves to the case q = 2; 
the dynamical variables are given by [z^ : : z^: z^] G P^: 



-I- 



(6a) 
(66) 

I • (6c) 

(6d) 

Notice that the renormalization map is symmetric for the exchange of the special 
state with the other state. In figure [71 we provide a plot of the spontaneous 
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Figure 7. Ising model on the spider web: spontaneous magnetization vs. z 
variable for z G M for various values of interaction. Note that the transition is 
rather gentle; this is due to the fact that the density of zeros is low near the 
transition point 



magnetization vs. interaction that confirms the presence of the paramagnetic phase for 
all antiferromagnetic interactions and of the ferromagnetic phase for all ferromagnetic 
interactions. 

5. Sierpinski gasket 

We can generate the Sierpinski gasket by infinite iterations of the decoration shown 
in figure [HI 




Figure 8. Some iterations of the decoration which generates the Sierpinski gasket. 

The dynamical and physical spaces are the same as in section 21 in this case each 
side of each triangle is counted just once, so we have no multiplicity issues. We record 
for sake of completeness the renormalization map for any value of q: 

= zl • ((/ - 3) - 2) ((/ - 1) + 3 • • z| ■{q-2){q-l) + 

4-{q-l) + 3-4-z^-{q-2){q-l)+ (7a) 

+5-4-{q-l) + 5-z^-4-{q-l) + z^ 
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iTp = ^1 ■{q-4)iq- 3) (9 - 2) + (z| + 2 • 2p • ^|) (q - 3) {q - 2) + 

+ (2 • 2p • z| + ■ 2|) (g - 3) (g - 2) + 

+3 • • 2| • (9 - 3) (9 - 2) + (2p • z| + 2 • 4 • ((/ - 2) + 

+ (2 ■ ■ ^1 + 2 • 4 • -^1 + 2 ■ 4) (9 - 2) + 

+ (2™ • + 2 • 4 • ^i) (9 - 2) + 
+ (2 ■ 4 • 2| + ) (9 - 2) + 

+ (2 ■ • 2p • 2| + 4) (9 - 2) + 4 • (9 - 2) + 



(76) 



2,2 



(7c) 



^ = zl-{q-5){q-4){q-3) + 3 (z| + 2 • Zp ■ z|) (g - 4) (g - 3) + 
+3 • • z| • (g - 4) (9 - 3) + 

+3 (^2| + 2 • 2p • z| + 3 • 4 • -^l) (9 - 3) + 
+3 (zp ■ z| + 2 • 4 • Z|) (q - 3) + 
+3 (2™ • z| + 2 • 4 • Z|) - 3) + 4 • (9 - 3) + 

+8 

Once more, if g' = 2, the third equation decouples and again we obtain a map defined 
on a 

Xm = 4:- z^- z^- z^-^ z^ (8a) 

^ = • (3 • + Znn) . {Sb) 

This is equation 3.2 in [11] or equation 14 in [12]. The exact and numerical results 
(figure [9|) tell us that we have no phase transitions at finite temperature; we have 
zeros in the thermodynamical limit only for T = in the physical domain. The 
paramagnetic fixed point is attracting for all points in the positive real axis, so that 
we cannot have a ferromagnetic phase. This behaviour is similar to that of the linear 
chain. Again it is easy but tedious to compute the map for the general case in presence 
of an external magnetic field; we will give here the exact expression for the Ising case: 





- 


3 • hX • • 4) + • 4 








(9a) 






• • Z^ -\- h^ha • {z^ -\-2 ' Z^ 


Zrm 












h^h^ • • Zrm -\- 2 • z^- z^^ -\- 


hi- 


Zrm 


4 


(96) 






■zL-ZBn + ^>hu ■ (4 + 2 • ^ff 


Zrm 




+ 








hX • (4 • ^- + 2 ■ 4 • ^fc) + 


hi- 


Zrm 




(9c) 


6^ 


+ 


+ Khl ■ {z^ ■ 4) + 
{zff ■ 4) + /^^ 4 








{M) 
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Figure 9. Fisher zeros for the Sierpinski gasket. The left column shows zero es 
in the interaction coordinates, while the right column show zeros for the inverse 
interaction coordinates; the three rows correspond to different q = 2,3, 4. 

Notice that, once again, we have complete symmetry for exchange of the special state 
with the other. Since we have no ferromagnetic phase, the Lee- Yang zeros do not 
accumulate to the positive real axis, as shown in figure [TOl 

6. Cylinders 

In this final section we provide an example of non-uniform lattices, i.e. lattices in which 
several types of edges are used. We present a lattice obtained as the quotient of the 
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Figure 10. Standard field coordinates: Lee- Yang zeros for the Sierpinski gasket. 
The picture is for z = 1.3 but quahtatively depicts all ferromagnetic interactions. 
As we expected, zeros do not accumulate on the real positive axis, and their 
structure is quite complicated as it is made of pieces with genuine zeros and 
pieces with points of the indeterminacy set. 

square lattice with a translation. Such lattices can be regarded as being generated 
by decorations in figure [11] and figure [121 For these lattices, a very special case of 
non- uniform lattices, we recover results that can be found in a completely equivalent 
way using the transfer matrix method; in this framework the transfer matrix is indeed 
the renormalization map. 

In each one of these lattices we have two type of edges. One is a regular 2-edge and 




Figure 11. Decorations that generate the skewed cylinders Z^/r • Z(l, 1) for 
r = 3,4,6, respectively. The cylinders are obtained by substituting infinitely 
many times the triangle (square, hexagon) with the corresponding decoration 




Figure 12. Decorations that generate the cylinders Z^/r • Z(l, 0) for r = 3, 4, 6, 
respectively. The cylinders are obtained by substituting infinitely many times the 
triangle (square, hexagon) with the corresponding decoration 

the other is respectively a 3,4 or 6-edge. Consider, for instance, the simplest lattice 
in figure [m i.e. the 3-skewed cylinder and let q = 2. The dynamical variables are the 
Boltzmann weights [z^ : Zp], [z^ : z^] G x P-^. The renormalization map will leave 
pair interactions (i.e. the second factor) invariant and will induce a 3-spin interaction 
on the first factor according to the following formula: 

^nm = Zaz2 ' Z^ -\- 3 • Z^ (^Z^ • -\- Z^ • ) + 

+ z^ ■ z^ (10a) 
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^ = Zan- Z^- Zg + 3 ■ Zgi (z^ ■ Z^ + Z^ ■ Zg) + 

+ z^-zl- Z^ (106) 

ir„ =z^ (loc) 

iTg =^e (lOd) 

Notice that we can arrange the map as a hnear map in the order 3 variables, parametric 
in the order 2 variables: 

^^\_f z^ + z^ 3{zt-z^ + zl-z^) 



[ ^S' )~[ {zi-z^ + zl 



zg, 3{z^-z^ + z^-z^) / V 



whose corresponding matrix is the transfer matrix of the system. Since we are deahng 
with a projective space we can factor out the polynomial 2|^-2;g +4- 2;g {ifz^+ZQ ^ 0), 
and defining: 



a{[z^ : Zq]) 



zt + zt 



■ [z^-zi + zi-4y 

we can rewrite the matrix in the much simpler form: 

a 3 
1 3 

Computing a in the standard interaction coordinates C (for C ^ ±i), we obtain: 



c= - 







^ 1 




H ' 1 ' 1 ' 

\ 


' 1 ' 1 ' V- 

- / 







Figure 13. Standard interaction coordinates: The set of non-analyticity points 
of the Green function; it looks hke a subset of the Fisher zero set of the I? lattice. 

We can compute the Green function in the variable ^ and obtain the set depicted 
in figure [13] for the non-analyticity locus. 

Notice that in this case the matrix of degrees does not satisfy the Perron- Frobenius hy- 
pothesis of |Appendix A.2[ in fact, the matrix is parabolic, i.e. it is not diagonalizable, 
with generalized eigenvalue 1. Therefore, we have to use a variation of the argument 
that we provided; the convergence of the Green function in this case is much slower 
(logarithmic) and the plot looks less definite. Indeed, one can obtain the set in an 
analytical way; in fact, one can easily check that in this case the appropriate version 
of the Green function is proportional to the logarithm of the norm of the maximal 
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eigenvalue of the matrix. The non-analyticity locus is therefore contained in the set 
of points such that we have at least two eigenvalues with maximum norm. Such a 
condition is easily expressed in an analytic form and the resulting set agrees with the 
numerical picture. 
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Appendix A. Numerical study of rational maps 
Appendix A.l. Fixed points and basins of attraction 

In cases where it is easy to locate all stable fixed points of the map (e.g. for maps on the 
Riemann sphere C), it is possible to obtain their basins of attraction in the following 
way. First of all we find a stable neighbourhood of each fixed point i.e. a ball of small 
radius such that its image is contained in itself. Having fixed a maximum number 
of iterations, we discretize a bounded region of the physical space in a finite number 
of pixels^ and for each pixel, we apply iteratively the map starting from the center 
of the pixel until we end up in one of the aforementioned stable neighbourhoods. If 
this happens in the given maximum number of iterations, we declare the initial pixel 
to belong to the attracting basin of the corresponding fixed point and we color it 
according to the speed of convergence (the faster the lighter) and to which fixed point 
it is attracted to. If the point does not fall in any neighbourhood in the given number 
of iterations, it is coloured black. 

Pictures obtained in such a way show the unstable set of the map as the boundary of 
the basins of attraction; moreover, showing which points are attracted to which fixed 
points, the pictures provide basic information on the asymptotic dynamics of the RG 
map. 

Appendix A. 2. Green Function 

We recall that a rational map on a multiprojective space ^ lifts to a polynomial map 
that is separately homogeneous in each factor, i.e.: 
/ : P^i X • . • X P^p ^ P^i X • • • X P^p 

{[z(% • • • , ^ ([F(i) (zW, • • • , , • • • , [F(^) (z(i), • • • , .(^))]) 

where each i^^*^ is such that: 

j 

and d'j is the degree of F*^*^ with respect to z'^^^ Considering d'j as an integer- valued 
matrix we can find its eigenvalues; in good cases we expect (via the Perron- 
Frobenius theorem) a simple real maximal eigenvalue p+ > 1 such that its associated 
(normalized) eigenvector W-^ has all non- negative coordinates. In such cases we can 
define the Green function: 

^ = lim — - (^+,log 

n^oo p'j^ \ 
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Notice the similarity of this function with the free energy of the system. In fact, for 
hierarchical lattices we have the following expression for the free energy: 

^= lim / ^ log ^o(f^(z^'\'-'.z^pA) , 
n^oo #edges Tn V V / / 

where is the partition function of the starting hypergraph. Moreover, if we call 
Si the number of edges of type i that belong to the starting hypergraph, we can 
express the total number of edges of the n-th approximation to the hierarchical lattice 
as ^i j Si {D^Yj. For generic (5, this expression is obviously asymptotic to p^. As 
explained in [13j there are results that state that the two functions ^ and ^ are equal 
in the uniform case with mild assumptions on but there is no general result for 
the no n- uniform case. Moreover notice that in the uniform case the matrix D is just 
a number, therefore most of the computations are made easier. 

We remark that we can exploit the homogeneous nature of the map to obtain a clever 
(and geometrically converging) way of numerically computing the Green function. In 
fact, let us define the sequence Zn of normalized iterates and the sequence An of the 
corresponding norms as follows: 



a; 



■0 — 
we can write: 



^0 



/K 



3 

^MI(^0 



Therefore, iterating the previous expression we get: 

n-l 



n n 



Taking the logarithm and considering logA^ as components of a vector logA^ in a 
p-dimensional space and again as elements of the p x p- matrix we obtain the 
following expression: 



log||(F-)({z^^)}) 



= logAn +I^logAn-l +1^ logAn-2 



-D" logAo. 



When we compute the scalar product with the maximal eigenvector of D we are 
projecting on the corresponding eigenspace, therefore the expression can be rewritten 
as: 



(^+, log An) + (^+,L>logAn-l) + 
= (^+, log An) +P+ (^+,logAn-l) 



(^+,DMogAo) = 
'^pI (^+,logAo). 



Dividing by the normalization term we get the following expression for the Green 
function: 
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that is geometrically convergent (if p+ > 1) and can be computed numerically with 
very good approximation as the are bounded. As a last remark notice that in the 
uniform case the expression reduces to 



where d is the degree of the map and A^s are just numbers. 
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